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Using post-Newtonian equations of motion for fluid bodies that include radiation-reaction terms 
at 2.5 and 3.5 post-Newtonian (PN) order (0[(v/c) 5 ] and 0[(v/c) 7 ] beyond Newtonian order), we 
derive the equations of motion for binary systems with spinning bodies. In particular we determine 
the effects of radiation-reaction coupled to spin-orbit effects on the two-body equations of motion, 
and on the evolution of the spins. For a suitable definition of spin, we reproduce the standard 
equations of motion and spin-precession at the first post-Newtonian order. At 3.5PN order, we 
determine the spin-orbit induced reaction effects on the orbital motion, but we find that radiation 
damping has no effect on either the magnitude or the direction of the spins. Using the equations 
of motion, we find that the loss of total energy and total angular momentum induced by spin-orbit 
effects precisely balances the radiative flux of those quantities calculated by Kidder et al. The 
equations of motion may be useful for evolving inspiraling orbits of compact spinning binaries. 



I. INTRODUCTION AND SUMMARY 
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. This is the third in a series of papers on motion and gravitational radiation in the post-Newtonian approximation 
to general relativity (GR). In Paper I 0, we developed a method of direct integration of the relaxed Einstein 
equations, whereby we wrote the Einstein equations as an inhomogeneous fiat-spacetime wave equation together with 
ly-j ' a harmonic gauge condition and solved them formally in terms of retarded integrals over the past null cone of a given 
field point. We obtained formal solutions both for the far-zone gravitational waves and for the near-zone fields needed 
to obtain equations of motion. We then developed the near-zone fields in a post-Newtonian expansion, in powers of 
q-i e ~ (v/c) 2 ~ Gm/rc 2 , where v, m and r represent typical velocities, masses and separations in the system, and G and 
c are the gravitational constant and speed of light (both set equal to unity henceforth). Each power of e represents 
^Jy one "post-Newtonian" (PN) order in the series (e 1 / 2 represents one-half, or 0.5PN orders). The near-zone metric was 
evaluated through 3.5PN order in terms of instantaneous, Poisson-like integrals over distributions of perfect fluid. 

In Paper II 0] we specialized the equations of motion to binary systems of non-rotating, suitably spherical balls 
of pressureless fluid, whose size is small compared to their separation, and derived the two-body equations of motion 
through 2.5PN order and including the 3.5PN terms (calculation of the 3PN terms is ongoing). Through 2PN 



order, the equations of motion are conservative, and our results were in complete agreement with calculations of 
others^, 0, la, IE The 2.5PN and 3.5PN terms represent the effects of gravitational radiation damping and 
their post-Newtonian corrections. We showed that they lead to losses of orbital energy and angular momentum that 
correctly match the losses calculated from the gravitational-wave flux to infinity U\m They are also consistent with 
a recent calculation of 3.5PN terms in the equations of motion using the post-Minkowskian approach 11]. 

In addition to the formal question of deriving suitably well-defined equations of motion in general relativity, a 
problem that dates back to the earliest days of the subject, this work is motivated by practical considerations. The 
operation of a network of ground-based, kilometer-scale gravitational-wave interferometers, and the research and 
development toward a space-based interferometer have made it critical to obtain highly accurate theoretical models 
for the evolution of and gravitational wave emission from two-body systems. One of the leading candidate sources 
for detection by interferometers is the inspiral of a binary system containing black holes or neutron stars, and the 
preferred method of detection, optimal matched filtering, requires theoretical "template" waveforms that are accurate 
to fractions of a wave cycle over the potentially thousands of cycles in a "chirp" signal whose frequency and amplitude 
increase with time as the binary inspirals. This is why calculations to high PN order are needed. 
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However, the two-body equations of motion of Paper II (and indeed of all high-PN order calculations done to 
date) have a shortcoming - they assume non-spinning bodies. But spin may play a critical role in binary inspiral, 
particularly involving black holes. Spin-orbit and spin-spin coupling leads to precessions of the spins of the bodies and 
of the orbital plane, the latter effect resulting in modulations of the amplitude of the gravitational waveform received 
at a detector. Furthermore, spin effects contribute directly to the gravitational waveform, and to the overall emission 
of energy and angular momentum from the system [12i Il3| . 

From the point of view of gravitational-wave data analysis, spin complicates matters by increasing the number of 
parameters that must be estimated in a binary-inspiral matched filter; this can significa ntly decrease the accuracy 
with which any parameters can be measured, and can increase the computational burden (14, Ha . Ha. flTl . Il8| . 

Thus, it is desirable to have as complete a theoretical picture as is reasonable of the effects of spin in relativistic 
binaries. Formally, spin effects first enter at the 1PN level. By analogy with quantum mechanics, the spin orbit 
contribution to the energy is of order L • S/r 3 , where L is the orbital angular momentum. This can be re-expressed 
roughly in the form 5E ~ (mrv)(mRV)/r 3 ~ (m 2 /r)(R/r)vV ~ e£jv, where R and V are the radius and rotational 
velocity of the spinning body, and En ~ m 2 /r is the Newtonian orbital energy Indeed, the 1PN effects of spin 
have been derived by numerous authors from a variety of points of view, ranging from formal developments of the GR 
equations of motion in multipole expansions |2Ct |21| , to post-Newtonian calculations 22] , to treatments of linearized 
GR as a spin- two quantum theory, with the concomitant spinning-body interaction potentials p3l l24j . For a review 
of these various approaches, see |25| . 

The resulting two-body equations of motion can be written in the form 



r n + a PN + a S o 



(1.1) 



where a = ai — &2 is the relative acceleration, and where the 1PN point mass and spin-orbit contributions are given 
by 
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where x = Xi - 
i = (m 2 /mi)Si 



X2, r = |x|, n = x/r, m = m\ + m^, rj = mim2/m 2 , v = vi — v 2 , r = dr/dt = n • v, S = Si + S2, 
f (mi/m2)52, and Ln = x x v. The 1PN equations of spin precession are given by 
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These equations assume a specific definition of the center of mass and spin of each spinning body, given by xi = 
m \ X Ji P* x d 3 x, and Si = Si(l + \v\ + 3m,2/r) — |vi x (vi x Si), with analogous expressions for body 2, where p* 
is the so-called "conserved" or baryonic density [Eq. I|2.1() ]. mi is the baryonic mass, and Si = J x p*x x vd 3 x, is the 
baryonic spin, with x = x — xi . The spin orbit terms in Eq. H.2b ]l and the spin precession equations (|1.3[) derived 
by our method are in complete agreement with earlier derivations [25| . 

Turning now to radiation-reaction effects, it is clear that no spin contributions can occur at the 2.5PN order 
corresponding to the leading-order damping effects. The simplest way to see this is to recall that, in a certain gauge, 
the 2.5PN radiation reaction acceleration may be expressed in the form 



i 2.5PN 
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5 dt 5 



(1.4) 



where i lJ is the trace-free moment of inertia of the system, and where we sum over repeated spatial indices. But 
since, to lowest order, f ^ contains no velocities, it cannot contain spin terms at lowest order, even if the bodies are 
rotating; in other words, spin terms will be of higher, 3.5PN order. Similarly, in calculating Si from x x a2.5PN it is 
also clear that no spin terms can arise at lowest order. Hence, the leading contributions of spin in radiation reaction 
must be at 3.5PN order. 

We therefore make use of the 2.5PN and 3.5PN equations of motion for fluid systems derived in Paper I and 
calculate the equations of motion and spin precession for two spinning bodies. These are our central results: the 
2.5PN contributions to radiation reaction are given by the standard "point-mass" result 

(1.5) 
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and the spin-orbit contributions to radiation reaction are given by 
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The point-mass, 3.5PN contributions are as given in Paper II, Eqs. (1.2) and (1.3d). We also find the radiation-reaction 
contribution to spin precession 
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(1.7) 



Although the latter result arises effectively from a cancellation of many terms, it should really come as no surprise. 
Our bodies are assumed to be axially symmetric, and as such, should not couple to gravitational radiation on the basis 
of their rotation alone. Their precessions could lead to gravitational radiation and to back-reaction if, for example, 
they were rotationally flattened, but those precessions and the flattening are themselves lPN-order effects, and so 
should result in radiation reaction effects of at least 4.5PN order, beyond the level at which we are working. 

We have verified that these equations lead to a loss of total energy and total angular momentum that matches 
precisely the energy and angular momentum flux at infinity from binary systems of spinning bodies, as calculated by 
Kidder aZ. |I2,E3l- 

In the conventional terminology of spinning bodies in general relativity, the center of mass of each body is specified 
by a so-called "spin supplementary condition" ( Appendix lAl. given by a parameter fcssc whose value is typically 
zero, one or 1/2. Our center-of-mass definition corresponds to the value fcssc = 1/2- Appendix iBl displays the final 
equations of motion corresponding to the "covariant" SSC given by fcssc = 1- 

The remainder of this paper provides details. In Sec. [H]we derive the equations of motion and spin precession to 
1PN and 2.5 PN order, and demonstrate explicitly the absence of observable radiation-reaction spin effects at 2.5PN 
order. In Sec. IIIII we complete the derivation to 3.5PN order, and verify the agreement with energy and angular 
momentum flux. A variety of technical matters and detailed equations are relegated to Appendices. 



II. POST-NEWTONIAN AND 2.5PN EQUATIONS OF MOTION 

A. Foundations 



We wish to analyse a binary system consisting of balls of perfect fluid that are sufficiently small compared to 
their separation that tidal interactions (and their relativistic generalizations) can be ignored, but that are sufficiently 
extended that they can support a finite rotational angular momentum, or spin. At Newtonian order, the result 
is essentially trivial: the equation of motion for body 1 is d 2 x.i/dt 2 — — ni2x/r 3 + 0(mR 2 /r 4 ), where R is the 
characteristic size of the bodies. Spin plays no role whatsoever, because the Newtonian interaction does not depend 
on velocity. But at post-Newtonian order, there are velocity-dependent accelerations of the schematic form mv 2 /r 2 , 
and thus, taking into account the finite size of the body and expanding about its center of mass, one could anticipate 
acceleration terms of the form (mVR)v/r 3 ~ Sv/r 3 . However, the combination of finite size and spin introduces an 
ambiguity in the definition of the center of mass of each body. At Newtonian order, the center of mass is defined 
naturally by = f. pxcPx, where p is the density and where the integral is over the body in question. But 
at post- Newtonian order, depending on the choice of "density" used, there could be correction terms in the center of 
mass of the form mT 1 j pv 2 xd 3 x ~ vS/m. Because "center of mass" is an arbitarify chosen point, these differences 
have no physical content in the end, as long as results are expressed in terms of measurable quantities (such as the 
total energy of the system, or the radar range to the surface of the body), but they do result in equations of motion 
with different explicit forms. This has given rise to the concept of "spin supplementary condition" (SSC), a statement 
about which center of mass definition is being used; this concept is discussed in Appendix [S] 
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Because we intend to work to post-Newtonian orders that include 1PN, 2.5PN and 3.5PN terms, we need to define 
centers of mass and spins in a way that is physically reasonable, calculationally simple, and easily transformable 
to definitions based on other criteria. In post-Newtonian theory, the simplest "density" available is the so-called 
"conserved" density, given by 

p*=p^u°, (2.1) 

where p is the mass energy density as measured by an observer in a local inertial frame momentarily at rest with 
respect to the fluid, g is the determinant of the metric, and u° is the time component of the fluid four-velocity. 
Making the reasonable assumption that p is directly proportional to the baryon number density, then, by virtue of 
conservation of baryons, we have that (pu a )- a = 0, where semicolon denotes covariant derivative, and the index a is 
summed over all spacetime values. This is equivalent to the exact continuity equation 

d P */dt + V • 0*v) = , (2.2) 

where v % = u % /u° is the ordinary (coordinate) velocity of the fluid. The use of p* in defining such quantities as center 
of mass is convenient because of the fact, based on the equation of continuity, that 



d_ 

dt 



J p* (t, x')/(x, x')dV = J p* (t, x')v' • V'/(x, x')dV , (2.3) 



for any suitably regular function /(x, x') and for integration over a complete body. 

Accordingly, we will define the baryonic mass, center of mass and baryonic spin of each body in our system to be 



A 
-1 



m A = I p*d 3 x, (2.4a) 
/ p*xd 3 x, (2.4b) 

J A 



S A = e zjk / p*x j v k d 3 x, (2.4c) 

J A 

where x J = x J — x A and v J = v J — v J A . For future use we will also define a tensorial spin quantity 

S l l =2 [ p*x [ ' l v ]] d 3 3 



5 A —l \ p x l v J 'a x . 
'a 



SJ = e^ k S k A , S A = ^ k S A k , (2.5) 

where [] around indices denotes antisymmetrization. We will demonstrate in Appendix lAlthat this definition of center 
of mass corresponds to the SSC value fcssc = 1/2. With these definitions, the baryonic mass tua is constant, and the 
velocity, acceleration and rate of change of spin of body A are given by 

v A = m^ 1 / p*vd 3 x , (2.6a) 

J A 

sla = m^ 1 I p*ad 3 x, (2.6b) 

J A 

dS\/dt = e ljk [ p*x j a k d 3 x. (2.6c) 

J A 

Note that, because J A p*5td 3 x = 0, we do not need the "barred" version of the acceleration in Eq. 



B. Baryonic equations of motion and spin precession 



We begin by working to 1PN and 2.5PN order, reproducing a number of well-known 1PN formulae for spinning 
bodies, and establishing some results that will be useful when we go on to 3.5PN order. Since we are only interested 
in radiation-reaction aspects of spin, we can ignore the 2PN terms in the equations of motion; these produce only PN 
corrections to the spin equations of motion. Incorporating the 2PN terms would require a more accurate definition 
of spin, including PN corrections. Such a framework has not been developed to date and is beyond the scope of this 
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paper. Initially, we derive everything in terms of our baryonic definitions of center of mass and spin, then later, we 
transform to more relevant definitions. 

We use the continuum equations of motion derived in Paper II, Eqs. (2.23), (2.24a) and (2.24c), with all quantities 
expressed in terms of the conserved density p* . They are given by 



where 

a PN 
a 2.5PN 



cPx'/df = U' 1 + a PN + 4. S p N , (2.7) 



v 2 U' 1 - 4vV'[/ J - 4Z7Z7' i - 3i/Z7 + AV l + 8v j V [i ' j] 

+l* > i-*£ + \x' i , (2-8a) 

o (5) , (5) (4) (3) 4 (3) (3) 

\x 3 (T 3 --5 13 T kk ) + 2v ] I 13 +2V 3 T j +-U> Z l kk -X' l3k V k , (2.8b) 

DO O 



where commas denote partial derivatives, overdots denote partial time derivatives, and (n) above quantities denotes 
the number of total time derivatives. The potentials used here and elsewhere in the paper are given by the general 
definitions 

J l x - x l 
X(f) = J p*(i,x')/(i,x')|x-x'|dV, 

Y(f) = |p*(t,x')/(i,x')|x-x'|W, (2.9) 

with specific potentials given by 

U = E(l), V i = E(u') , 

$! = E(« 2 ) , <ff = S(«^ fe ) , $ 2 = E(C/) , 

X = X(l), X l = X(v l ), X 1= X(v 2 ), X 2 =X(U), 

Y = F(l). (2.10) 

The multipole moments of the system Z y , X yfe , and J qj , and additional moments that will be relevant at 3.5PN 
order are defined in Appendix |0 I n Eq. (|2.8b(l , we have eliminated terms that are purely functions of time whose 
main effect is to generate an overall center of mass motion for the system that is irrelevant for our purposes. This is 
discussed in Appendix iDl 

We now multiply the equation of motion l|2.7|l by p* and integrate over body 1, expressing the variables x and v as 
x = xa + x and v = va + v, where A = 1, 2, depending on the body in which the point lies. To get the acceleration 
of body 1, we divide the result by mi. We use Eqs. I|2.4I) and I2.6fl to simplify where possible. We expand the various 
potentials in powers of x/r, and keep terms proportional to v x x, as well as "internal" terms proportional to v 2 , 
and p* If. The random part of the "w 2 " terms can be thought of as pressure contributions to the internal structure, 
so we no longer are treating the bodies as purely pressureless. We make extensive use of virial relations derived in 
Appendix [E] to simplify expressions dependent on the internal structure of each body. We also assume that the bodies 
are approximately spherically symmetric in their own rest frames (we ignore centrifugal and Lorentz-boost flattening, 
which will be higher-PN-order effects), so that integrals of internal expressions with an odd number of spatial indices 
(corresponding to odd-Z spherical harmonics) can be set to zero. We also ignore internal body terms that vary as x 2 , 
which represent quadrupole and higher "tidal" effects, and which vanish as the bodies' sizes tend to zero. 

The Newtonian term gives = —m2x\ 2 /r 3 , where, in this paragraph, we denote x\ 2 = x\ — x\, r = |xi2 1 , and 
n = X12/?*. To illustrate the approach, we show two examples of the terms that arise at 1PN order, namely 



— [ p*v 2 U^d 3 x = — — { p*{v 2 1 +2v 1 -v + v 2 )d 3 x 
mi Jx mi Ji 



p*'{x — x'Y ^ 3 , m 2 x\ 2 m 2 xi{S l 3 — d'tf'n?) 
|x — x'| 3 r 3 r 3 



2v{ 3l 2 x\ 2 27j m 2 x\ 2 S[ K v k m 2 (^ - 3nV) 

/t-i — m 2 v l — ; , \ Z - LL S 



mi r' mi ^ mi 
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and 



/ p*UU' i d 3 x 

mi J x 



/ p*d A x 

mi J x 



3 „ m 2 m 2 _ 
-a x H ^x • X12 



m 2 x\ 2 



m 2 x^ (S^ — Ztfn 3 ) 



fii m 2 x\ 2 ^1 m 2 x\ 2 m 2 x\ 2 



mi 



mi 



(2.12) 



where Ti, ftf , ^i and TC X are defined in Appendix [E] These entirely internal, structure-dependent terms can be 
eliminated via virial relations; however they can and do generate spin terms at 1PN and 3.5PN order. Detailed 
expressions for these virial relations may be found in Appendix [E| 

In the combination of 1PN terms 4V l + in Eq. 1)2. 8a|) . the time derivatives generate accelerations inside the 
potentials. To the order needed for our purposes, we must therefore substitute the Newtonian and 2.5PN continuum 
terms for those accelerations and carry out the same procedures for the integrals as described above. Those Newtonian 
and 2.5PN acceleration terms inserted into 1PN expressions will then generate respective 1PN and 3.5PN point mass 
and spin-orbit contributions. 

Turning to the 2.5PN terms, Eq. 



(|2.8bjl . the integrations lead to no explicit internal or spin terms; however for 
our 3.5PN accurate work, the multipole moments and their time derivatives must be evaluated to 1PN order, and in 
those 1PN corrections, spin terms will appear (see Appendix . 
The resulting equation of motion for body 1 is 



m 2 7l 



+ (ai)pN + K)sO + K)2.5PN 



(2.13) 



where 



01) 



PN 
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5PN 
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m 2 
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2(v k 2 S[ k -v k 1 S 2 ,k )-l(v k 1 S{ k 



,fc djk\ 



m 2 



,3 



(<F fe - 3n j n k ){ Vl - v 2 y{-S\ k + 2S lk ) . 



q (5) i (5) (4) -, (3) (3) 

-x{(Ti --^ X kk ) + 2v{ V 1 --^-n 1 l kk -^n l v?n k V k 



3 r 2 



(2.14a) 

(2.14b) 
(2.14c) 



where Sj[ = Sj /mA, and where the "mass" m 2 in the Newtonian term is now given, to PN order, by the baryonic 
mass plus \$l 2 , which corresponds to the total mass-energy of the body. 

We calculate the precession of the spin in a similar manner. Starting with dS\/dt = e yfe J 1 p*x : >a k d 3 x, we expand 
about the baryonic centers of mass, keeping terms that depend on x x v, dropping internal terms with odd numbers 
of spatial indices, and throwing away terms that vary as R 2 or higher. Notice that, even in Newtonian theory, such 
(m 2 /r 3 )(miRl) terms occur, and represent standard quadrupole coupling; in the Earth-Moon system, these lead to 
the precession of the equinoxes. However, as we wish to deal with compact bodies, we shall ignore such effects. At 
1PN order, the only terms in Eq. I|2.8a|l that contribute are those that have explicit velocity dependence. The result, 
at 1PN order is 



r<?M - 7712 



55j(vi • n) - 35j(v 2 • n) - 2n 4 (Si • v) + (« x - 2v 2 ) l (S 1 ■ n) 



(2.15) 



Strangely, however, there is a 2.5PN contribution to the spin equation of motion, arising from the u-dependent term 
in Eq. I|2.8b|) . and given by 



(4) (4) 

(^iksPN = si i" si t> 



(2.16) 



However, we notice that, since the spin is constant to lowest order, then, to 2.5PN order, the right-hand-side of Eq. 
(|2.16|) is a total time derivative, and thus can be moved to the left-hand-side and absorbed into a redefinition of the 
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spin. In any case, it is clear from the argument made in Sec. [I]that this is a spurious effect, because there is a gauge in 
which the 2.5PN terms in the equation of motion have no explicit velocity dependence. This 2.5PN gauge, sometimes 
called Burke-Thorne gauge 2&1 > ma y be obtained from our gauge by the coordinate transformation 

(3) i (3) 
X 1 -> X 1 - X 1 V 1 +-X 1 1" , 

3 

2 (2) 

* t + -T" . (2.17) 

Applying this transformation to our definition of baryonic spin Eq. I|2.4c(l . and recalling that p*d 3 x is an invariant 
quantity, we find that 

(3) (3) 

(^Burke-Thorne = S{ - S{ 1" +S{ 2* . (2.18) 

Hence, absorbing the 2.5PN terms into the definition of the spin is equivalent to defining our spin in the Burke-Thorne 
gauge. Although this redefinition eliminates 2.5PN terms from the spin precession, it will have consequences at 3.5PN 
order. 



C. Post-Newtonian conserved quantities 

It is useful to verify that, at 1PN order, the equations of motion and spin precession we have found admit suitable 
conserved quantities for energy, angular momentum and momentum or center-of-mass motion. Taking the equations 
of motion (|2.13() . i|2.14a|) and l|2.14bjl . contracting with miVi, summing over both bodies, and using the equations of 
motion to extract time derivatives (see Appendix [FJ, one obtains a conserved total energy. Doing the same procedure 
with a cross product with mixi, and combining with the spin precession equations, one obtains a conserved total 
angular momentum. Finally, multiplying by mi, summing over both bodies and extracting time derivatives, one 
obtains an expression for the system center of mass. These lPN-conserved quantities are given by 

1 2 1 17111712 3 4 3 2 m 2 1 ITllTlK 

E = - mi „ 1 --__ + - mi „ 1 + - mi „ 1 _ + -_ r - 

1 m 1 m 2 , . , . 

-(7v 1 -v 2 +v 1 -nv 2 -n) + (1^2), (2.19a) 



1 2 1 m 2 \ 1 TOi77l2 



4 r 

mi(xi x vi) ^1 + |ui - ) ~ ' ~'"x - v , + v.j • n x, - n) 

1 2 , r, m 2\ 1 



+Si ^1 + -vi + 3-iJ - - Vl X (VI X Si) 

+H!l(2n-\^) x(nxS x ) + (1^2), (2.19b) 
r \ 2 r J 

I = mixi ^1 + X -v\ - + l -v x x Si + (1 ^ 2) , (2.19c) 

where, to the 1PN order needed, the masses shown are the total masses of each body, given by vn\ + ifii and 
i7i2 + ^^2- The first thing to notice about these conserved quantities is the absence of a spin-orbit contribution to 
the total energy. This well-known result is merely a consequence of our choice of SSC. Converting from our SSC to 
the "covariant" SSC, for example, gives the standard spin-orbit term (see Appendix [Bj . These conserved quantities 
can also be derived from global definitions of energy, momentum and angular momentum, as discussed in Appendix 
mi It is straightforward to verify, by taking a time derivative of Eqs. I|2.19|l and substituting the 1PN equations of 
motion l|2.14aj) . (|2.14b|) and spin precession 1)2.15(1 . that E, J and I are constant to 1PN order. 



D. The proper spin 

The total angular momentum J, Eq. (|2.19b|l . has been written in a form that appears to have an orbital piece, 
plus 1PN corrections, a spin piece, plus 1PN corrections, and a final, spin-orbit piece. Although the split is somewhat 
arbitrary, it is useful in that, if we identify the "proper" spin of each body by the collection of "spin" terms in Eq. 
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(|2.19b|) . then Eq. (|2.15() is equivalent to the standard spin-precession equation (|1.3|) . In addition, we found that, in 
our gauge, there was a 2.5PN contribution to the spin precession [Eq. 1(2.16(1 ]. but since that contribution was a total 
time derivative to 2.5PN order, it could be absorbed into a redefinition of the spin. We therefore define the proper 
spin of each body to be 

1 m,\ 1 ( 3 ) - ( 3 ) 



S\ = S{\1 + ~v( + 3^j - -[vi x ( Vl x Si)] 1 - S\ T» +S{ I» , 
S 2 = (1^2), (2.20) 
With this definition, the spins Sa satisfy Eq. 1(1. 3|) . with no 2.5PN contributions. 

E. System center of mass and transformation to relative coordinates 

Choosing our coordinates so that the "center of mass" quantity I vanishes, and defining the transformation from 
individual to relative coordinates to 1PN order by 

x\ — — -x % + Sx 1 , 

TO 
777 -1 

x\ = -x*+Sx\ (2.21) 

TO 

we obtain 

5 j = 1 S JH _ ™W _ i^-e«V(5* - e) , (2.22) 
2 to V r ) 2 dm 

where 5m = m\ — mi . These transformations do not affect the Newtonian term in the acceleration, of course. In the 
1PN and spin-orbit terms they will only produce 2PN effects, which we ignore. The multipolc moments that appear 
in the 2.5PN terms in the equation of motion 12.14c|) must also be converted to relative coordinates, keeping any PN 
and spin-orbit corrections generated by Eqs. 12.2111 : this is treated in Appendix IU1 In addition, in the 2.5PN terms, 
multiple time derivatives of the multipole moments will generate accelerations, for which the 1PN relative equations of 
motion including spin-orbit terms must be substituted; in explicitly 3.5PN terms, the Newtonian equation of motion 
suffices. 

Calculating the relative acceleration a = ai a.2 using Eqs. 1)2.13(1 and 1)2. 14(1 . and converting to relative coordinates, 
we obtain the equations of motion given by Eqs. (|1 . 1|> . (|1.2(l and 1)1.5(1 . 

In terms of these relative coordinates, the energy and angular momentum of the system to 1PN order including 
spin terms then take the form, 

„ ( 1 TO 3 , „ a 1 , . n TO 1 TO „ 1 /TO\ 2 1 

E = M { ^ 2 ~ 7 + g(l - 3 ^ 4 + 2 ( 3 + ^ 7 + 2V + 2 (7) } ' (2 - 23a) 
1 / 2 ( x m "I t i 1 n 



J = /iL N il + -(1 - 3 V )v 2 + (3 + 77)- +S + -^nx[nx(45 + 3^)] , (2.23b) 
[2 r J 2 r 

where /x = r\m is the reduced mass. Notice that we do not keep the 2.5PN contribution to J arising from the 
conversion from baryonic spin to proper spin using Eq. 1(2. 20(1 . since E and J are only well defined up to 2PN order. 
In Sec. IIII CI we will use these expressions together with the 3.5PN equations of motion to compare E and J with the 
corresponding fluxes of radiation to infinity. 

III. 3.5PN EQUATIONS OF MOTION 
A. Equation of motion 

To obtain the 3.5PN contributions to the equations of motion including spin terms, we take the 3.5PN fluid 
expressions shown in Appendix [Dl multiply by p* , and integrate over body 1. We follow the same procedure as in 
Sec. Ill Bl expanding potentials about the baryonic centers of mass of the bodies, keeping point-mass terms, internal 
self-energy terms, and spinlike terms (terms linear in xv) and discarding tidal-like terms. Many terms in Eq. I(D4() 
make only point-mass contributions, such as terms of the form r 2 x k {d 7 I tk /dt 7 ), x^U{d 5 T^ /dt 5 ), v 2 v l (d 4 I kk /dt 4 ), 



9 



UU' 3 (d 3 T lJ /dt 3 ), and so on, because they do not have the proper mix of x and v to generate a spin. However terms 
such as r 2 v k (d 6 T k /dt 6 ), x 3 U{d 4 T^ /dt 4 ) or <S» l (d 3 l kk /dt 3 ) will generate spins. 

Some of these terms generate contributions of the form of an integral that is purely internal to body 1 multiplied by 
a multipole moment; examples include T^ 3 x k d 5 l jk /dt 5 , £l{ k v\d 4 l jk / dt 4 , Ti{ v k d 3 I jk / 'dt 3 . Virial relations must then 
be applied to these internal integrals, to sec if any spin terms arise. But in the virial relation l|E3|l involving T^ J and 
f?/, spin terms occur one PN order higher, hence there will be no contributions from these terms to the acceleration 
at 3.5PN order. On the other hand, the virial relation l|E4|) involving H. lJ does have a spin contribution at the same 
order, so those virial-induced terms must be kept. 

To lowest order, the mass multipole moments X 4 - 7 " and their derivatives do not contain spin terms, however the 
current moments J 13 " do contain both point- mass and spin terms at lowest order. 

In addition, the combination of 1PN terms 4V l + hX'* in Eq. (|2.8a|l . will generate accelerations whose 2.5PN terms 
will produce 3.5PN point-mass and spin terms that must be included. We must also re-express the 1PN spin-orbit 
terms of Eq. I|2.14bjl in terms of the proper spin of Eq. I|2.20[l : the 2.5PN contributions there will generate 3.5PN terms 
in the equation of motion. Finally, in the 2.5PN accelerations of Eq. (|2. 14c|l . we must include the 1PN corrections 
to the multipole moments as well as the 1PN terms in the equations of motion that are generated by the many time 
derivatives; these corrections will also contain spin contributions (Appendix Q . 

The result for the 3.5PN acceleration of body 1 is an expression too lengthy to reproduce here. The point mass 
terms reproduce Eq. (4.2) of Paper II, apart from small differences resulting from our use of the gauge transformation 
(|D2fl to remove purely time dependent terms from the 2.5PN and 3.5PN equations of motion. After transforming to 
relative coordinates and obtaining the relative acceleration, our point-mass terms match the corresponding expressions 
of Paper II precisely. The expression involving spins is equally lengthy. 

Calculating a\ — a\ using Eqs. (|2.14a|) . (|2.14b(l . I|2.14c|l and our 3.5PN terms, substituting Eqs. I|2.21|l and the 
time-derivatives of the multipole moments l|C2() . and expressing the spin results in terms of the total spin S and the 
spin quantity we obtain the final relative equation of motion terms as given in Eqs. (|1.2I) . 11.5fl and i|1.6[) . 



B. Spin precession 



We now want to calculate the precession of the proper spin S\ to 3.5PN order. A time derivative of Eq. i|2.20|) gives 



,m 2 



S\ 11 + vf + >; 
+S{ ( 2v x ■ ai - 3 



1 



r j 2 



r 



«i( v i-si) 

-la^vx-SO-^Kai-SO 



(4) (4) (3) (3) 

-S[ '£'■' +S{ T J -S{ 1" +S{ T j 



(3.1) 



We repeat the method of Sec. Ill HI to determine the contributions of 3.5PN fluid terms to the time derivative of 
the baryonic spin Si, by calculating e %] J 1 p*x 3 a k 5PN d 3 x. Only terms in a3 5PN [Eq. I|D4|) ] that have explicit v 
dependence will contribute a spin term. Notice that, as we have discussed, the 2.5PN contribution to S\ cancels the 
relevant terms in the last line of Eq. (|3.1|) . For ai, which appears in the 1PN terms in Eq. I|3.1|l . we must substitute 
the 2.5PN equations of motion; for S\ in the final 2.5PN terms in Eq. (|3.1|l we must substitute the 1PN precession 
equations; finally we must use Eq. (|2.20() to convert from Si back to the proper spin Si to the appropriate order. 

The result is the 1PN spin precession of Eq. (|1.3|) . plus a lengthy 3.5PN expression. However, using the fact that, to 
lowest order Si = 0, together with the identities listed in Appendix [0 it is straightforward to show that our lengthy 
3.5PN expression is in fact a total time derivative, given by 



(S 



1J3.5PN 



rfm d 
5 ~dt 

+nSi • n 



Si 



(14 - 13aV - (10 + 15a)r - -(27 + Ik*) — 

3 r 



15(4 - a)v 2 - 25(4 + 3a)r 2 - — (17+ 37a) — 

3 r 



-— nSi • v 
— -^vSi ■ n 



(18 - 14aV - 6(9 + 10a)f 2 - -(34 + 161a) — 

3 r 

(26 - 20a)v 2 - (66 + 45a)f 2 - ^(50 + 59a) — 

3 r 



- vSi • v(40 + 25a) 



(3.2) 
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where a = mi/ mi. As such, it can be moved to the left-hand-side and absorbed into a meaningless, 3.5PN term in 
the redefinition of the spin. As a result, we find, to little surprise, that radiation reaction makes no contribution to 
the precession of the spins [for a physical justification, see remarks following Eq. (f 1 . T|> ] . 



C. Comparison with fluxes of energy and angular momentum 



The fluxes of energy and angular momentum in gravitational waves from a binary with spin-orbit interactions were 
derived by Kidder et al. 0, 0| , and are given by 



dE _ . 
~~dt ~ 
dJ ■ 



- Ego 
Jso , 



(3.3) 



where we include only the lowest-order "Newtonian" and 1PN spin-orbit contributions. After transforming from 
Kidder's &ssc = 1 formulae, Eqs. (3.25a), (3.25c), (3.28a) and (3.28c) of Ref. to our /esse = 1/2 using the 

transformation x — > x — (v x £)/2m, we obtain for the fluxes, 



E N 
Eso 

Jn 
Jso 



8 rj 2 m 
15 r 4 
8 rj 2 m 3 



{12v 2 - llf 2 ), 



Ln • S 



[rir 2 



15 r 6 
-L N • £ (l8r 2 - 19v 2 
8 7y 2 m 3 



37, 2 - 12™) 
r / 



m 
I— 

r 



5 r 3 



2v 2 - 3f 2 - 2— 



4 r\ vn J 2 to 



5 r 3 

TO 

-I n x 

r 

+rv x 



3 r 



(^ 2 



1 TO 



)(<S - £) - o f — n x (7v x 5 + 5v x |) 



17 



(n x 5) ( 6f 2 - y v 2 + 2 



29 to 



3 r 
m 
r 



(nx5) h 24v 2 - 30f 2 |5(nx() h 4^ - 5 



3 r 



(n x |) ( 6f 2 - 6v 2 



5 TO 



4 TO 

37 



3 r 



+V X 



^LnLn 



, , n o 23?7l\ ^ f o 29 , TO 

vxS 18f 2 - I2v 2 + v x £[ 15r 2 v 2 - 7— 

3 r J \ 3 r 

9 92to\ f i o 59 to 

l 8w 2 _ _ + £ 20r 2 - 13w 2 - — - 
3 r J \ 3 r 



S 30r z 



(3.4a) 

(3.4b) 
(3.4c) 



(3.4d) 



We now calculate the time derivative of the energy and angular momentum expressions (|2.23f) . and substitute 
the equations of motion, including 1PN, spin-orbit, 2.5PN and 3.5PN spin orbit terms, along with the 1PN spin 
precession equations (recall there are no 3.5PN contributions to the spin precession). After recovering the fact that 
all 1PN point-mass and spin-orbit contributions cancel, leaving E and J conserved to that order, we find that the 
changes in E and J due to 2.5PN and 3.5PN spin-orbit radiation reaction are obtained from the following expressions, 



E — j«V • (a 2 .5PN + a3.5PN-so) 7 
J = px (a 2 . 5 pN + a 3 . 5PN _so) ■ 



(3.5) 



Initially, the results do not match the flux expressions above. However, by making use of the identities listed in 
Appendix [0 we can show that the difference between the expressions in all cases is a total time derivative. These 
can thus be absorbed into meaningless 2.5PN and 3.5PN corrections to the definition of total energy and angular 
momentum. Thus we have established a proper energy and angular momentum balance between the radiation flux 
and the evolution of the orbit, including spin-orbit effects. 



IV. CONCLUSIONS 



We have derived the equations of motion for binary systems of spinning bodies from first principles, including the 
effects of gravitational radiation reaction, and incorporating the contributions of spin-orbit coupling at 3.5PN order. 
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We found that the spins themselves are unaffected by radiation reaction. The resulting equations of motion are 
instantaneous, dynamical equations, and do not rely on assumptions of energy balance, or orbital averaging. They 
may be used to study the effects of spin on the inspiral of compact binaries numerically. We have focussed attention 
on effects linear in the spins, corresponding to spin-orbit coupling; the effects of spin-spin coupling can in principle 
also be calculated with our approach. These issues will be the subject of future work. 
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APPENDIX A: SPIN SUPPLEMENTARY CONDITIONS 



When we deal with systems containing spinning bodies, the fact that the bodies have finite size introduces an 
ambiguity in the definition of each body's center of mass. This has given rise to the concept of "spin supplementary 
condition" (SSC) [201 l2l| . which is a condition whose role is to fix the center of mass. One defines the antisymmetric 
tensor S 1 ^ , given by 

S^ = 2 [ (x^-x [ %)T^°d 3 x, (Al) 

J A 

where T a @ is the source stress-energy pseudo-tensor that appears in the "relaxed Einstein equations" given by Oh a P = 
-1671-1-^, and is a combination of the stress-energy tensor of matter and terms quadratic in the gravitational fields 
h al3 (see Sec. II. A. of Paper I), and where x\ is to be identified with the world line of the center of mass of the body. 
The center of mass is then fixed by imposing the following condition: 

Si°-fc ssc S>i = 0, (A2) 
where kgsc typically has the values 1, 1/2 or 0. The value &ssc = 1 corresponds to the so-called "covariant" SSC, 

SjfllAii = 0. 

In this paper, our baryonic center-of-mass definition x\ = J A p*x l d 3 x corresponds to /esse = 1/2, which can be 
seen as follows. From Eq. (|G3a|) . we note that, to the appropriate order, r 00 = p*(l + \v 2 — \U). Calculating S t0 
following the methods outlined in Sec. [HI and including the spin term generated by J A p*x l v 2 d 3 x, we find directly 

that Sf = \S l lv A , and thus that fc S sc = 1/2. 

One can also show that the relationship between the centers of mass for each value of fcssc is given by 

(^(^(^W + ^sM)^- (A3) 
rriA 

For further discussion of the role of the SSC in PN calculations, see [ill. I25I l^ . 



APPENDIX B: EQUATIONS OF MOTION WITH SPIN IN THE COVARIANT SSC 



We now transform our key equations to a form in which the centers of mass of the bodies are defined by the 
covariant spin supplementary condition, fcssc = 1- That transformation is given by 

4 1/2) = XI + TT^Vi x S x , (Bl) 

2rni 

where the variables on the right-hand-side are in terms of /esse — li and where Si is the baryonic spin. For the 
relative coordinate, the transformation, to 1PN order, is 

x (1/2) = x + -!-v x £, (B2) 
2m 
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where £ is again the baryonic spin variable. Converting from the baryonic spin to proper spin and keeping only the 
2.5PN correction from Eq. I|2.2L)[1 (the PN corrections will not be relevant for this purpose), we have 



i 1 (3) i (3) 

c (1/2) = x + — v x £ - — v x £ T kk +— v x (e P k )e k , 
2m 2m 2m 



(B3) 



where £ now represents the proper spin variable. We substitute this transformation into the Newtonian, PN, spin-orbit 
and 2.5PN terms in the equations of motion, Eqs. ((1. 1(1 . 11.2(1 and (jl.6|) . keeping only PN spin-orbit terms, 2.5PN 
terms, and 3.5PN spin-orbit terms. Where accelerations arise, eg. in transforming the velocities in the spin-orbit 
terms or in the 2.5PN terms, we must employ suitably accurate expressions, in order to generate all appropriate 3.5PN 
spin-orbit terms. The result is to change the form only of the 1PN spin-orbit and 3.5PN spin-orbit terms, which now 
become 



aso 

a 3.5PN-SO 



1 



6-L N • (S + £) - v x (45 + 3£) + 3rn x (25 + £) }, 
r 1 



fjm J rn 

"5^ 



(l20v 2 



280t 2 + 453- 



Lw • S 



120v 2 + 280f 2 + 458— L N • £ 



150r 2 



66- ) L N ■ £ 



--rv x S (48u 2 
3 



15f 2 + 364- 



-n X S I 31t> 4 - 260?;-/ 



2^2 , 245f 4 _ 



irv x £ (29lv 2 
3 V 

689 o m 



705f 2 - 772- 



o m 4 m 2 
+ 537f 2 — + - — 
r 5 r 



1 



, A 00 / 869 o m m 44 to 2 

-nxC 115i> 4 - 1130u 2 f 2 + 1295f 4 v 2 — + 849f 2 — + - 

3 r r 3 r 2 



(B4a) 



(B4b) 



The spin precession equations are not affected by this SSC transformation, since we are working to linear order in 
spins. The energy and angular momentum in this SSC take the form 



„ [ 1 9 m 

E = u\-v 2 

^{2 r 

3, a 1, , o Tn 1 m o 

+ - 1 - 3r7 + - 3 + ?7 K- + -r/-r 2 
8 z r 2 r 



/iL N (l + i(l-3r/)u 2 + (3 + r7)- 
1 /j r m 



1 /m\ 2 



iii (™n x [n x (45 + 2£)] . - v X [v X £]) 
2m I r J 



(B5a) 



(B5b) 
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The spin-orbit contributions to the energy and angular momentum fluxes, from Eqs. (3.25c) and (3.28c) of 01 are 
given by 



Eso 



'SO 



r] 2 m 3 



15 r b 



Ln • S 



2,lv 2 - 12^ 



£ 51r 2 - 43u 2 



4 rfm 2 



5 r 3 

m 
— n x 

r 



- — (v 2 - r 2 )(S - £) - -r—n x (7v x S + 5v x £) 
3 r 3 r 



(n x 5) 



\ 3 r 



(nxfl (9f 2 -8« 2 --- 
3 r 



+rv x 

+v x 
1 



29 m 

(nx5)( y7 



2Av 2 



30r 2 + 5(n x £) -— + 4« 



5 m 



v x 5 18r 



-LnL 



S 30r 2 - 18u 2 



2 23 m 

" T7 

92 m 

T7 



3 r 



12w z -— I + v x £ I 18r 2 - — v 2 - 9— 

'' ' 3 r 



„2 5/ .2 



n n 71 m 

£ 35r 2 - 19w 2 

3 r 



(B6a) 



(B6b) 



By following the same steps as in Sec. IIII CI calculating the time derivative of E and J, substituting the relevant 
contributions from the equations of motion and spin precession, and extracting total time derivatives from the result 
(using kggc = 1 versions of the equations in Appendix IF1) . we verify that the losses of orbital energy and angular 
momentum are completely equivalent to the fluxes above. 



APPENDIX C: MULTIPOLE MOMENTS 



The multipole moments that appear in the radiation-reaction terms in Eqs. I|2.8bjl and l|D4|) are defined by the 
general expressions QG1[1 . Because the quadrupole moment X 1 ^ appears in the 2.5PN terms, it will be needed to 1PN 
order, while the remaining moments will be needed to only the lowest, "Newtonian" order. Substituting the relevant 
expressions for r"' 3 from Eqs. |G[3j!, and carrying out the standard procedure as in Appendix IU1 we obtain 



JiJ 


= m x x l l f 1 + ^ 


J-ijk... 


+ (1^2), 
= miXi "' + (1 


J 13 


= e io6 Lu^ j 


jijk 


= e iab (mxvlx^ 


M ljkl 


ij kl 1 

= miv{x 1 - -- 
1 



-4<\sf)+(1^2), 

+ |^" fc S 1 a)b )+(1^2). 

r 



+2x[ k S l ^ i vi ) + (1^2), (CI) 

where parentheses around indices denote symmetrization. Note that higher-order moments, such as Jii kl or M. • 7 - ys 
appear only in 3.5PN terms that were transformed away in Appendix^] so their explicit forms are not needed. 
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Converting to relative coordinates, using the 1PN correct transformation in the leading term of T %] \ we obtain 

V j = i)mi y (l + ~(1 - 3?7)v 2 - i(l - 2??)yJ + ^rjm 2 rS ij + r]x { ^ j)k v k , 
= -r)5mx ijk , 

v ikX = nm(i - iTi)x ljkl , 

= -^mZj^' - r 1 e iab x^A a '> b , 

= ^m(i - 3*7)Lj,a^ fc + -r}e iab x^> k ^ h , 

M ijkl = rpn{\ - 377) (v ij - x kl 

-^i 1 rn 2 r{n tj 5 kl + n kl 8 i3 - n^ k 8 1 ^ - n J W + 25 iik S^ j - 2S i] S kl ) 

+2r ? x (fe ^ )(l ^' ) , (C2) 

where Ln = x x v is the orbital angular momentum per unit mass, £ y = (ms/m^S^ + (mxjm^S^ and A u = 
m(Si j /m 2 -S\ j /mi) = (m/Sm)(e j - S ij ). 

Time derivatives of the moments may be calculated using the relative equations of motion in place of x l ; 1PN 
equations including spin contributions must be used in the leading term in X lJ , while Newtonian equations are 
sufficient for the remaining terms. Since spin effects in the moments and their time derivatives are already at 1PN 
order, the spins themselves may be treated as constants. 

APPENDIX D: 3.5PN TERMS IN THE FLUID EQUATIONS OF MOTION 

We start with the 2.5PN and 3.5PN accelerations terms in the fluid equations of motion, written in terms of 
the conserved density p* , Eqs. (2.24c) and (2.24d) of Paper II. We then make a coordinate transformation at 
2.5PN and 3.5PN order to eliminate the purely time dependent terms. These terms cancel when we compute the 
relative acceleration, but contribute to the gravitational-radiation induced recoil of the system. In other words, we 
choose coordinates that are fixed with respect to the recoiling center of mass of the system. This will simplify the 
transformation between the coordinates Xi and X2 of the individual bodies and the relative coordinate x, Eq. I|2.21|l , 
and will simplify our analysis of angular momentum (which can depend on the choice of center of mass) The required 
transformation is 



& X 2.5PN + ^ X 3.5PN 1 (Dl) 



where 



2 (3) 2 (2) 

S4.5PN = -T^--e^J«, (D2a) 

no (5) 9 (4) 1 (3) 

dXs5PN - 4200 1 75 6 J m M ■ [mb) 

However, because of the velocity dependence in the PN terms in the equations of motion, the change in velocity at 
2.5PN order, 8v l 2 5PN — bx\ 5PN will induce additional 3.5PN contributions in the equations of motion given by 

5al 5PN = 6v 3 2 5PN (2vW' 1 + v l U' 3 - US lj + V 3 ' 1 + X' ij ) 

+ \H.,pn(X^ -SUn. (D3) 

When we add these changes to the 3.5PN terms in Eq. (2.4d) of Paper II, we obtain the final 3.5PN contributions to 
the fluid equations of motion, given by 

1 (7) 

4.5PN = ^( 13r2xkSiJ ~ ~ X t X : >X k ) V k 

1 (6) 

+ (lOrVtf^ + 4{v ■ x)x k 5 t] - rV<P' fe - Ax t x 3 v k + 3x j x k v i ) V k 

oU 
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i (5) 

+ — (10(« • x)v k 6 ij - x k v 2 5 ij + 2(v ■ x)v l S 3k + 2x i v 2 6 jk - hx l v 3 v k + 4x j v i v k ) V k 
1 o 

1 (5) i (5) 

+ — (5r 2 [/ J - 35x 3 U + 59X^) T 3 + — {br 2 U l + 30x l U - 34X< 1 - 6x j X' ij ) I kk 

1 (5) 

+ (ISajJa;*^ - 6x 3 X' ik - 30x'X^ k - 15r 2 X> ijk - AY^ k + 5x l Y< i3kl ) V k 

1 (4) i (4) 

+ v *( v *S' k - v 3 v k ) l jk --(2x j U + 22V 1 - UX' 3 - \2v k X- 3k ) T 3 

+—{2Ax 3 v k U' i + \2v l x 3 U' k + \2x 3 V k ' 1 + bAv l X^ k - 72v 3 X' ik 
18 

(4) i (4) 

+12x 3 X' lk + QQX^ kl - 72X l ' 3k - hY' l3k ) V k (6v l U - 8V l - 2X- r ) I kk 

3 

(3) 

+ {2v 2 U- 3 - 8UU' 3 - 8v k V k > 3 + 3$f - 2<S>$ + X ' 3 ) V 3 

i 0) 

+ -(2v 2 U- 1 - 8v l v 3 U' 3 - 2AUU' 1 - 6v % U + 8V' 1 + \&v 3 V [l ' 3] + - + X' 1 ) l kk 



(48v 3 V k ' t - 12v i v k U' i - 6v 2 X' ijk + 2Av i v l X' 3kl - 24X i > 3k - 9X[ ijk + 6X? jk 



2 



(3) 

+\8v i X' 3k - 48v l X^ jk + 24UX> ijk + 24U' i X> jk - 18$i*'* + 6S' i (X' ik ) - Y> ijk ) l jk 

1 (7) 

+ (10x i x 3 8 kl - 25x k x l S ij - 9r 2 5 ij 5 kl ) V kl 

630 v ; 

i (6) 

+ — (4cM<f w - 2(v ■ x)5 t3 S kl - 10x k v l 6 ij ) l jkl 

1 (5) i (5) 

( v 2 S i3 6 kl + 6v l v 3 S kl + 5v k v l 5 i3 ) l 3kl UU5 3k + 10x 3 U k - hX' 3k ) l i3k 

45 45 

i (5) , (5) 

+ — (7X' i3 - 10x 3 U^) V kk + — (6x 3 X< lkl - Y' i3kl ) I 3kl 
45 54 

4 • (4) 
+—(U5 ij - v l U' 3 - 2v 3 U< 1 - V 3 -' 1 - X' i3 ) T 3kk 
45 



1 (6) Ifi (5) 

+ — (3r 2 e qtk + 2x l x 3 e q3k + 4x 3 x k e gij ) J qk -—x 3 v k e qjk J ql 
45 45 
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2 (5) 
+ — (2(v ■ x)t qlk - 2x i v j e qjk + 5x j v i e qik + 12x j v k e qij + 4x k v j e qij ) J qk 
45 

1 (4) 

(2Ue qik + 2x j U' i e qjk - Ax j U> k e qij + X' ij e qjk - AX^ k e qlJ - Ax 1 X< tjk e qlj ) J qk 

9 

o (4) A (4) 2 (3) 

+-(4vVe TO - v 2 e qlk ) J qk —x j U' k e qjk +-Ue qtk J qk 

2 (3) -I (7) 

— (vijj'i + 2v j U' i + V J > 1 + X^)e qik J qk x i l" kk 

9 840 

i (7) 1 (6) 1 (5) 1 (6) 1 (6) 
+ — x 3 JVkk + i jjjkk + _JJ,i JJjkk j qij jqkk J qik jqjk 

35 40 24 30 15 

i (5) (5) (5) 1 (5) 1 (5) 

_| _ v 3 ^jk jqik _ € qik jqjk _ e qij j-qkk^ ^_ x i J^kkjj ^^j" j^kkij 

15 30 15 

1 (4) 2 (4) 1 (3) 2 (3) 1 (3) 

— v i M kk " +-v j M ijkk +-U> 1 M jjkk +-U> 3 M i3kk — X^ k M ]kU , (D4) 
6 3 6 3 3 

where the relevant potentials are given in Eqs. I|2.10|l and the multipole moments are given in Appendix IU1 

APPENDIX E: VIRIAL RELATIONS 

Virial relations are statements about the internal structure of bound bodies that may be used to simplify the 
equations of motion. They assume that the body is either stationary or periodic over a suitable timescale. The 
simplest such virial relation, used in classical mechanics, states that 1/2 = 2T + H, where I = J p*\x\ 2 d 3 x is the scalar 
moment of inertia of the body, T is the internal kinetic energy, and f2 is the internal gravitational potential energy. 
For a stationary body, or when averaged over several internal timescales, we can set I — 0, and hence 2T + Q = 0. 
Here we generalize this to a variety of tensorial virial relations, including spin effects and effects at 2.5PN order. 

We define the tensorial moment of inertia of body A by 

ijf ee f p*x l xH z x. (El) 

J A 

Taking a time derivative and setting it to zero yields 1% — 2 f. p*x^v^d 3 x = 0. But from this and Eq. (|2.5|l . we can 
conclude that 

x l v 3 d 3 x = ^S l i . (E2) 

A second time derivative gives I 1 / = 2 f A p*v l v : >d 3 x + 2 J A p *x^ l a^d 3 x. For a J , we substitute only the Newtonian, 
post-Newtonian and 2.5PN terms from Eqs. I|2.7|) and 1)2.8(1 . In the equations of motion, these virial relations will 
be needed only to simplify PN terms, so for this application, only the Newtonian and 2.5PN terms will be needed. 
On the other hand, in the expression E = J M r m d 3 x for the conserved total energy to PN order, the virial relation 
will be needed to simplify the Newtonian contribution, and thus the 1PN spin contributions to the virial relation will 
be needed. Expanding potentials about the center of mass of each body as described in Sec. (|IIBJ) . keeping PN and 
2.5PN spin terms and 2.5PN "internal" terms, discarding terms that vanish as R 2 and higher as the body's size tends 
to zero, and setting time derivatives of the moment of inertia tensor to zero, we obtain the virial relation for body 1: 



2T\ 



f j + ^ = ^ (2n fc 5 1 fc(t (^ 1 ) -v£) + n^S{ )k {v k - 2v k )) 



, (4) , (3) (3) 

+S k{ * V )k -=-Sl\ j l kk -3Q\ jkl I kl , (E3) 
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The relevant internal quantities used here are defined below. A third derivative of f£ gives I a= 2 f. p*(3v^a^ 
^)d 3 x. The same procedure yields, for body 1, 



2H? j) - 



ij _ 3 ™2 q k(i j)k 
1 _ oTX* 5 ! n 



2 r 

1 (5) c (3) 

5 r d 2 

, (4) i (4) o (4) 

1 6 1 2 1 

_^( iaC W) w + 6JCf ij - 15/CI tfW ) 2*' , (E4) 



where we have used the lowest-order virial relation from Eqs. JE3|) and (E3J|, 2T[ J + Q{ = 0, and 2H\ J > - §/q 
— I^Sj n 3 ')* to simplify some of the 2.5PN order terms in Eq. I]E4[I. The relevant quantities used in these virial 
relations are defined by 

77 = ±J^p*v^d 3 x, 



, > >■ r— - w 



n y« _ 1 / / . (s - x' fjx - xjHx - xjHx - x'Y 3 3 
/ / p V ' V ' - (X = X ^ ( '~ y (g " X ' )3 d 3 xd*x> , 

^ s /■ [ ff> A*-W*-*) k i*-*Y Ma/ 

JaJa |x-x'| 5 



A J A 



s / / p y v ' ■ (* - - ^ - y d - ^) fc (^ - *')' Mx > , (E5 ) 



APPENDIX F: EXTRACTING TOTAL TIME DERIVATIVES 



Using the Newtonian equations of motion plus the 1PN spin-orbit terms, it is straightforward to establish a number 
of identities, which may be used to extract time derivatives from 2.5PN and 3.5PN terms. For any non-negative 
integers s, p and g, we obtain 



d 
di 



r q 



V 



2s-2 -p-1 



d fv 2a f p . 



dt V ri 



r q+l 

r q+l 
9 . 

v r 



pv 
pv 



4 - (p + q)v 2 r 2 2sr 2 ™ -pv 2 - + ^L N • (45 + 3®}, 
r r A r 



■ pv 
r r 



4 (p + q)v 2 r 2 - 2sr 2 - - pv 2 - + ■ (AS + 3£) 

2 r 6 



v rn 

2 



x (45 + 3£ 



(Fl) 
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Another set of identities, to be used only in 3.5PN terms, require only the Newtonian equations of motion: 

,,2s-2 -p-1 



d fv 2s rP . . 



m m 
pV 



....,= n I Ipv 4 — (p + q)v 2 r 2 — 2sr 2 

dt \ ri ) ri +1 ll v ; r r 

+2w 2 frx ( V ) j 

d fv 2s fP , A v ^-2 fP -i 



x x J 



r r 

i — \ \P vi — {p + q)v 2 f 2 — 2sr 2 



m 



(ft\ r« / ri+ l I r v ^ ~ r 



r 



2 • 

-2m— 



dt\r<> / r<? 



■P 1 f r 

i 1 pw 4 - (p + g)w 2 r 2 - 2sf 2 pv 2 — 

1 I L r r 

(t)V-\w)}. (F2) 



+u 2 fr 



APPENDIX G: TOTAL ENERGY AND ANGULAR MOMENTUM 

In this appendix, we express the global definitions of mass, momentum, angular momentum and center of mass in 
a PN expansion, including spin contributions. In Paper I, we defined the "source moments" 

P" = / T^d 3 x, (Gla) 
J M 



f T 00 x^d 3 x, (Gib) 

Jm 

f x aQ T 0b d 3 x, (Glc) 
Jm 

M ijQ = ( T ij x Q d 3 x, (Gld) 



I M 



JM 

where the capitalized superscript Q denotes a multi-index {xP = x 11 x 12 . . . x lq ) . The integrals are to be taken over 
a constant-time hypersurface M, which lies within the near zone, i.e. within one gravitational wavelength of the 
sources. The quantities P° and P l are naturally interpreted as the total mass-energy and momentum of the system, 
respectively, P = T l {q= 1) is the system dipole moment, which relates the center of mass to the origin of coordinates, 
and J 1 = J 1 (q = 0) is the total angular momentum. These quantities are conserved up to changes induced by a flux 
of gravitational radiation to infinity, that is, 

P" = - / T^d 2 Sj , 
JdM 

r = -e tab I T ih x a d 2 s ] , 

JdM 

r = P l - I T^x l d 2 S 3 . (G2) 

JdM 

where dM denotes the boundary of A4. Thus, these quantities are conserved and well-defined only up to and including 
2PN order. 

We now wish to evaluate these conserved quantities explicitly to 1PN order for a many-body system. Explicit 
expressions for r"^ through the relevant order may be obtained from Eqs. (5.7) and (5.9a) of Paper I, together with 
Eqs. (2.21) of Paper II which serve to convert expressions to dependence on p*. They are given by 

r 00 = p*(l + \v 2 -\u+ V + v 2 U + \u 2 + i$! - * 2 - 2v j V j + ]v j X'A , (G3a) 
y 2 2 8 2 2 4 / 



p*v j ^1 + ^v 2 -uj- ^X> j + ^-V fc ^8UV [k ^ ] + 4UU6 jk - X^U> k) + ^ k X< m U' m ^J , (G3b) 
p*vV + -L f / - U^WU 2 ) , (G3c) 



4tt V 2 
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where we have extracted total divergences of functions, such as V • (WVU) from the expression for t 00 , since they 
will not affect the integral quantities P° and P, to the orders considered. 

Expanding potentials about the baryonic center of mass of each body as described in Sec. (|IIBJl . keeping spin terms 
and other "internal" terms, but discarding terms that vanish as R 2 and higher as the body's size tends to zero, and 
applying the virial relations of Appendix lEl wc obtain the same expressions (I2.19|) for E, J, and I, as we found from 
the equations of motion. 
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